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Estimates  are  made  of  the  energy  losses  from  the  hot  channel  produced  by  (laser^ided)  elee> 
trie  discharges  in  the  atmosphere.  The  loss  processes  considered  are  heat  conduction,  coVtinuum 
'adiatio^,  and  line  radiation,^  For  the  conditions  calculated  to  exist  in  such  channels  20,000 
K,  2  X  10^^  enf’^),  Jw  total  enerv  lorn  is  found  to  be  in«U  compared  to  the  nraHni 
ener^  input  (r3  u/cm  along  the  length  of  the  channel).  Therefore  the  *stumpt!on  of  ne^Ugible 
losses  used  to  ^a'Ailate  the  channel  conditions  is  valid  and  the  expan  lion  of  the  hot.  channel  wiQ 
be  adiabatic  except  where  compressional  shock  waves  are  involved..  . 
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ENERGY  LOSSES  FROM  A  HEATED  AIR  COLUMN 


L  INTRODUCTION 

Experiments  at  NRL  which  use  the  reduced-density  channels  produced  by  laser-guided,  electric 
discharges  in  the  atmosphere  have  been  described.’  Initially  these  channels  are  narrow,  hot,  at  full 
atmospheric  density,  and  highly  overpressured.  They  expand  in  <  30  /us  to  form  atmospheric  pressure, 
quasi-stable,  reduced  density  channels.  Simulations  of  the  ohmic  heating  process  using  an  air  chemistry 
code-  have  assumed  radiation  transfer  and  shock  heating  may  be  neglected.  Hydrocode  simulations  of 
the  evolution  of  the  channel  densi'y  profile  have  also  assumed  the  expansion  process  taking  place 
behind  the  outgoing  shock  wave  is  adiabatic.^ 

We  estimate  the  radiation  loss  rate  from  the  channel  at  peak  current  using  the  air  chemistry  code 
predictions  of  the  channel  conditions.  We  also  estimate  absorption  lengths  both  inside  and  outside  the 
channel.  The  radiation  loss  mechanisms  which  have  been  considered  are  continuum  emission  due  to 
electron-ion  bremsstrahlung,  continuum  emission  due  to  electron-neutral  bremsstrahlung,  continuum 
emission  due  to  free-bound  transitions,  and  line  radiation.  Heat  conduction  losses  are  estimated  for  30 
fj.s,  i.e.,  during  the  whole  expansion  phase,  and  proven  negligible.  These  estimates  provide  reasonable 
support  for  the  assumptions  made  in  the  codes. 

II.  CHANNEL  CONDITIONS 

The  experir.iental  current  and  voltage  signals  for  the  electric  discharge  which  heats  the  channels 
were  simulated  by  combining  the  CH.VIAIR  air  chemistry  code^  with  a  simple  hydrodynamic  code  and 
the  external  circuit  equations.  The  predicted  expansion  was  consistent  with  Schlieren  and  image  con¬ 
verter  camera  photographs  of  the  channel  provided  suitable  initial  conditions  were  chosen.  Justifying 
these  initial  conditions,  and  showing  that  they  are  consistent  with  our  measurements  of  the  breakdown 
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process,  is  the  subject  of  ongoing  research  which  in  part  depends  on  understanding  the  radiati  e 
processes  occurring  in  and  around  the  channel. 

'  simulation  code  followed  electron,  gas,  and  vibrational  temperatures  and  various  species 
densities  as  well  as  the  radial  expansion.  The  values  taken  as  representative  for  the  first  current  peak 
and  used  to  estimate  radiation  losses  are 

T,  =  20.000  K, 

N,  =  1.6  X  10‘*  cm-\, 

-  6  X  10‘*  cm-^, 
radius  =  .25  cm 

where  N„  is  the  neutral  density.  The  dominance  of  Nil  and  Oil  lines  in  the  visible  spectrum  (Fig.  1) 
supports  these  estimates  of  the  conditions. 

During  the  next  29  /as  the  hot  gas  within  the  channel  expands  to  reach  pressure  equilibrium  with 
the  atmosphere  around  it.  At  equilibrium  the  measured  channel  radius  is  1.1  cm.  and  we  estimate' 
that  the  equilibrium  temperature  is  f,,  ~  Tj,  —  5000  K.  To  determine  an  upper  limit  i.f  the  conductive 
cooling  we  assume  a  time  average  channel  temperate  of  ~7000  K  and  a  radius  of  ~1  cm. 

in.  CONDUCTIVE  COOLING 

Determining  the  conductive  cooling  of  a  hot  body  generally  requires  solving  the  heat  conduction 
equation'* 

•|^  =  «V-r  (1) 

dr 

subject  to  initial  and  boundary  conditions.  In  this  equation  a  is  the  thermal  diffusivity  and  is  related  to 
the  thermal  conductivity,  k,  the  specific  heat,  c,  and  the  density,  p,  of  the  conducting  medium  by  the 
relation 


We  assume  that  the  hot  air  in  the  channel  can  transport  heat  quicxly  compared  with  the  cold  outside 
air  We  also  assume  that  the  heat  loss  from  the  channc*  will  prove  to  be  small.  We  may  then  consider 
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Fig  1  -  riie  iime  intriraicd  visib!^  spefirunr.  of  the  channel 


the  channel  to  be  a  heat  source  maintained  at  a  surface  temperature  of  ~7000  K.  It  is  not  necessary  to 
solve  Eq.  (I)  in  a  cylindrical  geometry  because  the  heat  from  the  channel  will  penetrate  only  a  short 
distance  in  30  ns.  We  may  make  a  simple  one  dimensional  argument  for  the  heat  loss  rate  per  unit 
surface  area/  If  the  heat  flows  outward  from  the  channel  to  a  characteristic  depth  8.  there  exists  a 
mean  temperature  gradient  of  T/8  (T  is  the  channel  temperature  which  is  large  compared  to  the  am¬ 
bient  temperature).  There  exists  a  corresponding  outward  heat  flow  per  unit  surface  area  of  kT/8.  The 
section  of  boundary  layer  in  contact  with  the  unit  surface  area  has  a  mean  temperature  of  T/2  and  is 


therefore  absorbing  heat  at  a  rate  of  pc 


8  where  6  is  the  rate  of  increase  of  the  boundary  layer 


thickness  and  the  rate  of  increase  of  the  volume  of  the  boundary  layer  section.  The  characteristic  value 
of  8  is  8/ 1.  We  may  equate  the  rate  at  which  energy  leaves  the  channel  surface  and  the  ••ate  of  energy 
increase  in  the  boundary  layer 


kT 


pc 


r  I 

T  t‘ 


(3) 


The  thickness  of  the  boundary  layer  after  time  t  is  therefore  given  by 

8^  “  2  a  r,  .  (4) 

(We  note  that  8  in  Eq.  (4)  is  just  a  factor  of  >/2  larger  than  the  characteristic  width  of  the  boundary 

layer  given  by  the  exact  solution  of  the  diffusion  equation;  see  for  example,  Carslaw  and  Jaeger,  "Con¬ 
duction  of  Heat  in  Solids,"  Clarendon  Press,  O.xford,  1959  p.  51.]  Substituting  the  value  of  thermal 

dift'usivity  for  ambient  air,  we  obtain  the  boundary  layer  thickness  at  39  ms, 

8  =  4  X  10'^  m.  (5) 

The  mean  heat  flow  per  unit  surface  area  is  then 

~  ~  5  X  10'>  — (6) 

8  sec  m* 

This  flow  rate,  sustained  for  30  ms  over  the  surface  area  of  a  unit  length  of  channel,  gives  a  total  energy 
loss  of  —  9  X  !0"-  J/cm,  and  is  negligible  compared  to  the  ~  3  J/cm  of  energy  ih  a  unit  length  of 
channel. 

IV.  LOSSES  DUE  TO  ELECTRO.\-ION  BRE.MSSTRAHLUNG 

A  summary  of  the  treatment  presented  in  Zel’dovich  and  Raizer^  for  deriving  the  "integrated 
emission  coefficient"  for  electron-ion  bremsstrahlung  from  the  "effective  radiation"  due  to  Coulomb 


scattering  will  be  given  (and  will  serve  as  a  guide  for  the  electron-neutral  and  bound-free  calculations). 
The  effective  radiation  is  defined  as  the  total  energy  emitted  between  frequencies  v  and  v  +  dv  when  a 
unit  electron  flux  passes  by  one  ion  (e.g.  one  electron  of  specified  velocity  and  moving  parallel  to  the  c 
axis  traverses  each  square  cm  of  an  x-y  plane  containing  the  ion).  The  effective  radiation  is  given  fiy 

(7) 


dq„  =  ~T-i= — —  dv 


3>/3  v^ 

where  Z  is  the  ion  charge,  m  the  electron  mass,  and  v  the  electron  velocity.  This  approximation  holds 


for  electron  velocities  satisfying  the  quasi-classical  condition 


/iv 


2nZe' 


-  «  1 


(8) 


which  is  true  in  normal  density  air  up  to  temperatures  of  10^  K.  To  account  frx-  the  actual  number 
of  electrons  in  the  gas  we  assume  they  have  a  Maxwellian  velocity  distribution 


f<v)  rf^v 


m 


lirkT] 


i/2 

•  g-mwVlkT  ^Jy 


(9) 


and  that  it  is  isotropic.  Each  unit  volume  then  contains  NJ'M  d\  electrons  moving  in  any  given 
direction  with  velocity  v.  Thus  there  appears  to  be  a  uniform  electron  flux  of  vA/^/(v)  d^\  in  that 
direction.  The  total  flux  impinging  on  an  ion  is  found  by  summing  over  all  directions. 

J  v.V^  /(v)  d^v  -  v/V^  /(v)  <\iTV^dv 


all  anijics 


=  vAC,  Flv)  dv 


where  F(v)  is  the  Maxwellian  speed  distribution 


Fiv)  dv  =  4it 


ni 


3/2 


y2  ^~'nW*^2kT 


(10) 


(11) 


2nkT] 

.Vlulliplying  the  effective  radiation  by  this  electron  flux  and  by  the  number  density  of  ions  gives  the 
energy  emitted  per  second  per  unit  volume  between  frequencies  v  and  v  +  dv  due  to  electrons  with 
speeds  between  v  and  v  +  dv. 


J,.  ^dvdv=‘  V  F{v)  dv  dq^iv).  (12) 

The  "spectral  emission  coefficient."  J,.,  is  defined  as  the  energy  emitted  per  second  per  unit  volume 

between  frequencies  v  and  v  +  dv  and  is  obtained  by  summing  expression  (16)  over  all  electron  speeds 

sufficient  to  create  photons  of  energy  hv 


(13) 


J^dv  =  N'^  Ng  dq/yv)  dv. 


The  minimum  speed  is  determined  by 


1  ■> 

^  hp. 


Substituting  from  relations  (7),  (11)  and  (14)  and  integrating  yields 


iZiL  .2?  .  g-hviKT 

3  ikTm  mc^  - 


The  integrated  emission  coefficient,  is  defined  as  the  energy  emitted  per  second  per  unit  volume 
and  is  found  by  summing  the  spectral  emission  coefficient  over  all  frequencies. 

,  C°°  ,  j  327r  [inkT  '  erg 

Je,  “  “  “T"  "1 —  — rr  K  — \ — •  (16) 

*'0  3  3/m  mc^h  cm^  sec 

For  the  conditions  in  the  channel  the  integrated  emission  coefficient  equals  —5.1  x  lO"*  J/sec  cm\ 


V.  LOSSEsS  DUE  TO  ELECTRON-NEUTRAL  BREMSSTRAHLUNG 


The  integrated  emission  coefficient  for  flectronmeutral  bremsstrahlung  may  be  found  by  starting 
with  the  effective  radiation  for  electron-neu  ral  collisions^ 

8  e^vV,„ 

where  (r„,,  is  '.he  momentum  transfer  cross  section.  A  Maxwellian  distribution  may  again  be  assumed 
for  evaluating  the  eivctron  flux.  We  sum  the  radiation  resulting  from  all  electrons  of  speed  exceeding 
v^in  and  multiply  by  the  neutral  density  (collective  effects  are  negligible  for  our  conditions)  to  find  the 
spectral  emission  coeftkieni 

1.^1"'  f  „S) 

The  evaluation  of  the  mt ’gral  ts  accomp'ished  by  repeated  integration  by  parts 


J  .v'  e-'-dx  «  +  I  f  C. 
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The  spectral  emission  coefficient  therefore  becomes 


_  e^N.N,  [j_U  '  +  ^  +  i]  ,-Tfj, 

3  c^  TTw)  2  A:r|  A:r  2)  ^ 

Integration  over  all  frequencies  gives  the  integrated  emission  coefficient 


80  N„N,  , 


3  h  yfir  I  w 


sec  cm-’ 


U'-mg  a  representative  value’  of  (r„^  (— 10~'^cm’)  we  find  the  integrated  emission  coefficient  to  be 
~2.4  lU’  J/sec  cm^ 


VI.  LOSSES  DUE  TO  FREE-BOUND  TRANSITIONS 


We  will  estimate  the  radiation  losses  from  the  recombination  of  elections  and  ions  by  treating  the 
bound  s’ates  as  hydrogen  like.  The  energj  levels  are  then  given  by 


-  — 


27r’  me*  Z’ 


The  classical  formula  lor  .he  effective  radiation  due  to  electron-ion  bremsstrahlung  (7)  may  be 
considered  as  describing  a  transition  between  hyperbolic  electron  orbits.  The  electron  motion  in  the 
higher  bound  levels  may  also  be  treated  as  quasi-classical.  Extending  relation  (7)  to  describe  a  transi¬ 
tion  between  hyperbolic  and  eliptic  ornns,  leads  to  an  estimate  of  the  electron  capture  cross  section  into 
a  giver,  level.* 

=  128  ;r’Z^ 

3-/^ 

The  energy  .'mined  per  second,  pei  unit  volume  from  transitions  by  an  electron  of  a  particular  speed 


into  a  given  level  is  given  by 


Substitutii’,.;  or  'r,„  yields 


J\.n  ~  lV''<r„,)(v,V^  Fiv)  dv). 


128  N-N,  Z* 

— 7-7= - 5 — 

3V3  m  o*v  li^iF 
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The  total  energy  emitted  per  second,  per  unit  volume  is  found  by  summing  over  all  speeds  and  all  lev¬ 
els.  (Extending  (23)  beyond  quasi-classical  limits  can  be  justified  by  comparison  to  quantum  mechani¬ 
cal  calculations) 


1287r^A("iV,Z^e">  p-  f(v)^v  «  ^ 

Inter-particle  interactions  lower  the  continuum  threshold  le?ving  only  a  finite  number  of  distinguishable 
levels.  The  summation  therefore  remains  finite  and  we  take  it  to  be  of  order  unity.  Substituting  for 
f  (v)  and  performing  the  integration  yields  the  integrated  emission  .coefficient  for  free-bound  transitions 


■'/A 


512  /V*/V,  Z*  e'®  kT 

3v^ 


m 


2nkT 


3/2 


^'lonizaiion  srg 
■cT  sec  cm^ 


(27) 


Substitution  of  the  channel  values  gives  ~8.6  x  10^  J/sec  cm‘. 


VII.  LOSSES  DUE  TO  LINE  RADIATION 


Tubulated  transition  probabilities  and  multiplicities^  allow  us  to  calculate  the  line  radiation  for  the 
dominant  species  NI,  NIL  01,  and  OIL  Initially  we  assume  all  lines  are  optically  thin.  We  let  a 
represent  the  quantum  numbers  of  the  upper  level  of  a  line  and  a'  the  quantum  numbers  of  the  lower 
level.  Tiic  number  density  of  the  upper  level  is  N„,  the  spontaneous  transition  probability  and 
the  frequency  of  the  line  is  The  power  emitted  per  unit  volume  in  one  line  is  then  given  by'® 


erg 


cm-^  sec 


(28) 


The  channel  will  be  in  a  state  of  partial  thermal  equilibriui.*  and  the  densities  of  various  excited 
states  of  the  ionic  species  will  be  determined  by  the  total  ion  density  N*  and  the  electron  temperature 
T" 

*  e 


.Y. 


-£■  kr 
N*g  e  "  ' 


cm 


(29) 


ZJT,) 

where  is  the  multiplicity  of  the  level  a,  is  the  excitation  energy  of  this  level  and  Z^iT^)  is  the 
partition  function  of  the  ion  as  a  function  of  temperature  and  is  available  in  tabulated  form.'^  The 


power  emitted  in  a  particulai  ion  line  is  then 


8 


(30) 


=  ^ - A-hv„ 


*aa  * ' aa 


2+(r,) 

and  the  total  power  radiated  by  all  lines  from  this  ion  spe>,:es  is  given  by 


_M:_y  „  ,_ES_ 

"“““cm' sec 


(31) 


A  similar  equation  with  A'®  replacing  N'^  and  Z^iT^)  replacing  Z^.(T,  gives  the  total  power  radi¬ 
ated  in  ihe  lines  from  the  neutral  atom. 


The  bulk  oi"  the  emitted  radiation  is  in  the  ultraviolet.  Summing  over  the  resonance  series  for  MI 
and  on  gives 

?  —  5  X  10^  J/sec  cm' 

and  over  the  same  series  for  01  and  NI 

P  ~  1  X  10^  J/sec  cm'. 

Vm.  REABSORPTION  OF  RADIATION 

Outside  the  channel  the  cold  air  will  be  quite  transparent  to  the  visible  frequencies  but  essentially 
opaque  below  —1800  A  due  to  the  Schumann-Runge  band  system  of  the  oxygen  molecule,"  (/—  0.1 
mm  for  A  —  1500  A).  Thus  if  significant  amounts  of  radiation  are  emitted  from  the  hot  channel  they 
will  be  absorbed  in  the  cold  gas  around  the  channel,  so  broadening  the  hot  channei. 

Within  the  channel  the  reverse  processes  of  inverse  electron-ion  bremsstrahlung.  inverse 
e.ectron-neutral  bremsstrahlung,  and  photo  ionization  will  each  contribute  to  the  opacity  of  the  channel 
and  we  may  estimate  the  resulting  frequency  averaged  absorption  coefficient. 

In  an  infinite  medium  in  thermal  equilibrium,  there  exits  an  isotropic  total  incident  power  of 
4tr7^  where  cr  is  the  Stephan-Boltzman  constant.  In  order  to  balance  emission  and  absorption  within 
the  medium  there  must  exist  a  frequency  averaged  absorption  coefficient  given  by 

+  lb 

- *rr 
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The  corresponding  mean  free  path  /]  *  l/ATj  may  be  evaluated  using  our  previous  values  for  the 
integrated  omission  coefficients  and  the  channel  temperature.  The  resulting  value  of  /|  ==  4  cm  justifies 
the  assumption  that  the  channel  is  optically  thin  to  the  continuum. 

Stark  broadening  parameters'"*  indicate  that  the  ultraviolet  ion  lines  will  be  collisionly  broadened 
to  —.05  A  and  that  they  will  have  absorption  lengths  of  —  2  x  10“^  mm  within  the  channel.  The  visi¬ 
ble  lines  will  be  Stark  broadened  to  widths  of  —4  A  corresponding  to  absorption  lengths  of  —1  cm. 
Visible  emissions  will  therefore  escape  the  channel.  Experimental  measurements  of  the  visible  line 
widths  (Fig,  2)  are  in  reasonable  agreement,  indicating  that  the  electron  density  estimate  is  good. 


I _ I _ I _ I _ I _ i_ _ „.J 

3980  3990  4000  4010 

X(X) 

Fig  2  —  The  profile  of  ihe  Ml  line  at  3995  A 
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IX.  LOSSES  DUE  TO  OPTICALLY  THICK  LINE  RADIATION 


We  have  seen  that  the  u.v.  lines  will  be  optically  thick.  Their  profiles,  as  seen  at  the  channel  sur¬ 
face,  are  broadened  and  saturated  at  the  black  body  level.  The  power  (emitted  per  frequency  interval 
and  per  unit  surface  area  is‘^ 

5(t-,r)  {1  -  c'l  (32) 

where  t  is  the  optical  depth  of  the  emitting  medium  and  S(i/,T)  is  the  Planck  function  integrated  over 

27r  steradians’^ 

'  1~' 

1  .  (33) 

The  Planck  function  is  appro.ximately  constant  over  the  width  of  the  line  so  the  total  power  emit¬ 
ted  per  unit  area  in  the  line  is 

/  S.  (34) 

We  must  evaluate  the  half  width.  A*',  of  the  profile  represented  by  Eq.  (32).  The  half  intensity 
points  occur  where 

r(j/  ±  It’)  =■  ln2.  (35) 

The  optical  depth  is  a  function  of  the  thickness  of  the  plasma.  /,  the  density  of  atoms  in  the  lower 
state  of  the  line,  .V„-,  and  the  absorption  cross  section  per  atom  irif). 

r(j')  =  (36) 
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The  half  width,  At',  of  the  optically  thick  line  then  follows  from 


ail'  ±  At^}  =  -r; — :  (38) 

where  we  have  combined  Eqs.  (35),  (36)  and  (37),  We  determine  for  each  u.v.  line  by  using  Eq. 
(38),  taking  F  from  the  Stark  broadening  data,  and  setting  /equal  to  the  column  diameter  ( /  =  .5  cm). 
Using  these  values  of  At'  we  evaluate  /  for  each  line  via  Eq.  (34).  Summing  over  the  Nil  and  Oil  lines 

and  multiplying  by  the  surface  area  of  a  unit  length  of  channel  gives  the  total  emission  per  unit  length 

Nil,  Oil  /»  =  1.5  X  lOM/seccm. 

The  density  of  neutral  atomic  species  is  considerably  less  than  the  ion  density;  however,  because  the 
atomic  levels  are  broadened  more  than  the  ionic  levels,  the  atoms  make  a  larger  contribution  to  the 
emitted  power, 

NI,  01  P=5x  10^  J/sec  cm. 

X.  CONCLUSIONS 

For  the  conditions  calculated  as  representing  the  electric  discharge  heated  channel  at  the  first 
current  peak,  the  emissions  from  a  unit  length  of  channel  [r  .25  cm)  due  to  various  radiative 
mechanisms  are  as  follow.s; 


Electron  Ion  Brem.  1  x  10''  J/sec  cm. 

Electron  Neutral  Brem.  5  x  10'  J/sec  cm. 

Recomb.  Radiation  2  x  10^  J/sec  cm, 

Line  Radiation  (thick)  6.5  x  (O'*  J/sec  cm. 

These  losses  are  small  compared  to  the  rate  of  energy  deposition  within  the  channel  (~- 10*  J/sec  cm)  at 
that  time.  Heat  conduction  losses  have  also  been  seen  to  be  small.  Thus  the  assumption  of  adiabatic 
'xpansion  is  reasonable  for  an  approximate  treatment  of  the  hvdrodynamic  expansion  to  pressure 
equiiibiu’m.  However  of  the  four  different  p'oeesses  of  radiation  emission,  one,  radiation  in  optically 
thick  lines,  is  strongly  temperature  dependent.  Therefore  if,  for  instance,  the  same  electrical  energy 
were  deposited  in  a  channel  of  only  one  half  the  present  diameter  (0  25  cm  instead  of  0.5  cm)  so  that 
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the  temperature  rose  to  —  80,000  K,  the  energy  radiated  in  optically  thick  lines  would  rise  more 
than  X 100.  Clearly  such  radiation  output  would  exceed  the  energy  input  and  the  hot  channel  would  not 
be  able  to  sustain  itself.  At  the  same  time,  the  emitted  radiation  would  be  absorbed  in  the  cold  air 
around  the  original  channel  in  an  annulus  <  1  mm  thick.  Thus  the  channel  would  rapidly  be 
broadened  to  a  diameter  of  —0.5  cm  and  the  temperature  reduced  to  a  value  at  which  radiation  losses 
were  no  longer  significant.  In  this  way  line  radiation  effectively  limns  the  temperature  to  which  air  can 
be  heated  by  an  electric  discharge  to  values  of  <  20,000  K  since  the  electrical  conductivity'*  varies  only 
as 
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